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1. Introduction 

One of the most appcahng and natural mechanisms for generating the tiny neutrino masses 
experimentally allowed is the (type I) see-saw mechanism [1]. This simple extension of the 
Standard Model (SM) with two or three right-handed Majorana neutrinos provides also a 
very attractive origin for the observed Baryon Asymmetry of the Universe, via leptogenesis 
[2,3]. A lepton asymmetry is dynamically generated in the out of equilibrium decay of the 
heavy Majorana neutrinos, and then partially converted into a baryon asymmetry due to 
(B + L)-violating non-perturbativc sphaleron interactions [4]. 

Unfortunately, a direct test of the see-saw mechanism is in general not possible. If the 
Dirac neutrino masses are similar to the other SM fermion masses, as naturally expected, 
the Majorana neutrino masses turn out to be of order 10^ — 10^^ GeV, so the right-handed 
(RH) neutrinos can not be produced in the LHC or future colliders, neither lead to other 
observable effects such as lepton flavour violating processes. On the other hand, although 
the see-saw mechanism can also work with Majorana masses as low as 100 GeV, which 
in principle are within the energy reach of the LHC, the smallness of the light neutrino 
masses generically implies in this case tiny neutrino Yukawa couplings and as consequence 
negligible mixing with the active neutrinos, leading also to unobservable effects. 

In order to obtain a large active-sterile neutrino mixing, some cancellations must occur 
in the light neutrino mass matrix. In this context, the small neutrino masses are not 
explained by a see-saw mechanism, but rather by an "inverse see-saw" [5], in which the 
global lepton number symmetry U{1)l is slightly broken by a small parameter, fi. In other 
words, the fine-tuned cancellations required in the light neutrino masses are not accidental. 
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but due to an approximate symmetry. The smallness of the parameter // is protected from 
radiative corrections (even without supersymmetry), since in the hmit /x — > a larger 
symmetry is realized [6] . 

Much attention has been devoted recently to this possibility, both in the context of 
LHC phenomenology [7] and in leptogenesis [8-10]. A consequence of the slightly broken 
U{1)l symmetry is the existence of two strongly degenerate RH neutrinos, which combine 
to form a quasi-Dirac fermion. This is interesting for leptogenesis, because it leads to an 
enhancement of the CP asymmetry, avoiding the strong bounds which apply to hierarchical 
RH neutrinos [11], and allowing for successful leptogenesis at much lower temperatures, 
T ~ 0{1 TeV) *. 

We focus here on a different scenario: we assume that the small lepton number violating 
effects responsible of light neutrino masses are negligible during the leptogenesis epoch, so 
B — L \s effectively conserved. This can be the case if global lepton number is broken 
spontaneously at a scale well below the electroweak phase transition [14] or, even if lepton 
number is broken at high scales, it leads to a CP asymmetry too small to account for 
the observed baryon asymmetry. The main difference with previous approaches is that in 
this framework the RH neutrinos combine exactly into Dirac fermions, and the total CP 
asymmetry vanishes. As a consequence, in order to generate a baryon asymmetry we have 
to rely on i) flavour effects and ii) sphaleron departure from thermal equilibrium during 
the leptogenesis epoch. 

Leptogenesis without neutrino Major ana masses, so-called "Dirac leptogenesis", has 
already been considered in the literature [15], but in a completely different set-up. In Dirac 
leptogenesis global lepton number remains exactly unbroken (except for the SM sphaleron 
interactions), so the light neutrinos are Dirac fermions made of the SU(2) doublet ul and 
the singlet RH neutrino vr. Realistic models contain not only the SM plus RH neutrinos, 
but also additional heavy particles to generate a non-zero lepton number for left-handed 
particles and and an equal and opposite lepton number for right-handed particles in their 
CP violating decay. 

This paper is organized as follows. In Sec. 2 we describe our leptogenesis framework, 
the CP asymmetries produced in the heavy Dirac neutrino decay and the basic requirements 
to generate the baryon asymmetry. In Sec. 3 we write the network of Boltzmann equations 
relevant for leptogenesis without Majorana masses. In Sec. 4 we present our results, both 
in the resonant and non-resonant regimes, and we conclude in Sec. 5. 

2. The Framework 

Our starting point is that above the electroweak phase transition the relevant particle 
contents for leptogenesis is that of the SM plus a number of SM singlet Dirac fermions Ni. 
Without loss of generality we can work on a basis in which Ni are mass eigenstates. In 
this basis the Lagrangian above the electroweak phase transition can be written as: 

C = CsM + iNj^Nj - M jNiNi - X^i P^da - X^'eaPRNji, (2.1) 

*See also [12, 13] for alternative (extended) models of leptogenesis with RH Majorana neutrinos at the 
TeV scale and without resonant enhancement. 
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where a,i are family indices (a = e,/i, r and i = 1,2,3,...), arc the leptonic SU{2) 
doublets, h = {h'^,h^)'^ is the Higgs field {h = iT2h*, with T2 Pauli's second matrix) and 
Pr^l are the chirality projectors. 

A quantitative illustration of the proposed scenario which can also account for the 
low energy neutrino phenomenology can be found in the context of the inverse see-saw 
mechanism. In this type of models [5], the lepton sector of the Standard Model is extended 
with two electroweak singlet two-component leptons per generation, i.e., 



4 = ] ,eRi,i^m,SLi ■ (2.2) 




In the original formulation, the singlets su were superstring inspired E(6) singlets, 
in contrast to the right-handed neutrinos ujn, which are in the spinorial representation. 
More recently this mechanism has also arisen in the context of left-right symmetry [16] and 
SO(IO) unified models [17]. 

At zero temperature the (9 x 9) mass matrix of the neutral lepton sector in the 
vl-i^'r-iSl basis is given by 

/ \ 



M = 



mjy M 



(2.3) 



\ M'^ fi J 

where mo and M are arbitrary 3x3 complex matrices in flavour space with 



mo = XaiV , (2.4) 

and V = 174 GeV being the Higgs vacuum expectation value. Moreover is a 3 x 3 complex 
symmetric matrix. 

The matrix M can be diagonalized by a unitary transformation, leading to nine mass 
eigenstates: three of them correspond to the observed light neutrinos, and the other six 
are heavy Major ana neutrinos. 

In this "inverse see-saw" scheme, assuming rriD,^ <^ M the effective Majorana mass 
matrix for the light neutrinos is approximately given by 

m„ = moM'^'^ iJ.M~^mj^ , (2.5) 

while the three pairs of two-component heavy neutrinos combine to form three quasi-Dirac 
fermions with masses of order M. The admixture among singlet and doublet SU (2) states 
(and the corresponding violation of unitarity in the light lepton sector) is of order mo/M 
and can be large [18-20]. This is so because although M is a large mass scale suppressing 
the light neutrino masses, it can be much smaller than in the standard type-I see-saw 
scenario since light neutrino masses in Eq. (2.5) are further suppressed by the small ratio 
n/M. Thus, in this scenario, for M as low as the electroweak scale the only bounds on 
the Yukawa couplings are those arising from constraints on violation of weak universality, 
lepton flavour violating processes and collider signatures [21]. 

It is important to notice that in the ^ limit a conserved total lepton number can 
be defined. This can be easily seen if, together with the standard lepton number Lsm = 1 
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for the SM leptons we assign a lepton number L^r = 1 to the singlets Ufn and s^j. With this 
assignment the mass matrix (2.3) with = conserves L = Lsm + Lj^. Then, the three 
light neutrinos are massless Weyl particles while the six heavy neutral leptons combine 
exactly into three Dirac fermions, Ni, which above the electroweak scale are given by: 



The smallness of the term can be easily understood if the total lepton number is 
spontaneously broken by a vacuum expectation value (a), with fi = f{a) [14]. In this case, 
light neutrino masses are a consequence of total lepton number being broken at an energy 
scale much lower than the electroweak scale {a) <^ v, and fi vanishes exactly at the heavy 
neutrino decay epoch. This scenario introduces one extra scalar singlet which couples with 
SL and the SM Higgs as 



and therefore can affect our results when considered in the framework of the inverse see- 
saw. In principle, there is a new Dirac neutrino decay channel, Ni Nja, which could be 
relevant for leptogenesis; in practice, as we will see the present mechanism only works for 
very degenerate heavy singlets. Mi ~ M2, therefore this channel is phase-space suppressed 
and our analysis will remain valid. 

In this framework, if jj, is effectively zero at the leptogenesis epoch, all processes con- 
serve B and L at the perturbative level. On the other hand, the sphalerons violate B+Lsm 
but conserve B — Lsm and, since the new heavy leptons are SM singlets, they do not change 
Ln. Therefore the SM sphaleron processes also conserve B — L. In brief B — Lis conserved 
by all the interactions of the model on scales above {a). 

Thus, one is effectively working in the limit in which the three light neutrinos are 
massless, the heavy ones combine into three Dirac fermions given by Eq. (2.6) and the 
relevant interactions are those in Eq. (2.1). 

Indeed, if leptogenesis occurs via the decay of heavy Standard Model singlets in a 
framework without Majorana masses above the electroweak scale, Eq. (2.1) has the relevant 
information. Thus our results will hold whether light neutrinos acquire masses by the 
inverse see-saw mechanism with the fi term generated at low scales as described above or 
by some other mechanism, as long as it does not imply the presence of new states relevant 
for leptogenesis. 

2.1 The CP Asymmetries 

One important ingredient that determines the baryon asymmetry generated in thermal 
leptogenesis in this scenario is the CP asymmetry produced in the decays of the heavy 
Dirac neutrinos iVj into leptons of flavour a, Caf- 



Ni = SLi + i^m . 



(2.6) 



(2.7) 



T{Ni ^ l^h) - V{Ni ^ IJi) 



(2.8) 



r(^i ^ ^ah) + T{Ni ^ IJi) ' 



a 
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(a) Tree (b) Wave function 

Figure I: The tree and one loop diagrams that contribute to the CP asymmetry in decays when 
the heavy neutrinos are of Dirac type. 



Because of the Dirac nature of N^, the only one- loop contribution to the CP asymmetry 
arises from the interference of the tree-level and the self-energy one-loop diagrams displayed 
in Fig. I, and it is given by [8, 22] ^ 

— 1 x-^ cii — 1 



^J2 (a- - 1T2 + (^^^)..V^V^ E VWiVWjP^p ' (2-9) 



where aj = M?/M? , gj = Tj/Mi and 



is the decay width of N^. In the last equation we have also introduced the effective mass mj. 
Notice that the above CP asymmetry is only part of the usual wave function contribution 
for Majorana Ni, so that in the Dirac case the total CP asymmetry exactly vanishes 



= Eem = 0, (2.11) 



a 



by CPT invariance and unitarity. Thus in order for leptogenesis to occur we must be in a 
temperature regime in which flavour effects [24, 25] are important. 

In writing Eq. (2.9) we have expressed the Yukawa couplings in terms of the projection 
coefficients Kai (related to the absolute values of the Yukawas) and some phases ^ai as 



Xai = y^iJ{\^X)iie'^-^ , (2.12) 



^We have already regularized the divergence at Mi = Mj using the resummation procedure of [8]. In [23] 

a different calculation was performed, obtaining that the regulator of the singularity is MjVj — MiVi instead 
of MiTj, but for the values of the widths that we consider (Fj 3> Fi), both results coincide. 
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where 




(2.13) 



Pal3 ~ Ppa ~ - 




(2.14) 



The factor ^^^'^\-^2^g-2 can be resonantly enhanced to Mj/2Tj if Mj — Mi ~ rj/2 [8]. In 
this regime the resonant contribution to the asymmetry becomes 



2.2 The Generation of B: Basic Requirements 

In the scenario that we are presenting B — Lsm — Ln is conserved by all processes. Since we 
want to explore the possibility of generating the cosmological baryon asymmetry exclusively 
during the production and decay of the heavy Dirac neutrinos, wc assume that at the 
beginning of the leptogenesis era B — Lsm = and Ln = (same abundance of and 
Ni)K Therefore it is clear that B — Lsm = after the heavy neutrinos have disappeared. 
If the sphalerons were still active after the decay epoch, the final baryon asymmetry, being 
proportional to B — Lsm, would be zero. 

Then, in order for leptogenesis to occur, one must be in the regime in which the 
sphalerons depart from equilibrium (which occurs at T ~ Tj) during the decay epoch. 
In this case, as described in next section, the baryon asymmetry freezes at a value Yb oc 
^b-Lsm(^ = Tj), which in general is not null. 

To go on we assume that the observed baryon asymmetry is generated during the decay 
epoch of the lightest heavy neutrino, Ni. For the sake of concrctcness we will assume that 
M3 ^ M2 > Ml so the maximum contribution to the CP asymmetry in Ni decays is due 
to Note that, in this scenario, for hierarchical masses the CP asymmetry is suppressed 
as (Mi/M2)^, instead of M1/M2 in the standard case with heavy Majorana neutrinos. This 
contributes to the impossibility of generating enough B with hierarchical heavy neutrinos. 
Moreover, if Ni and N2 have similar masses, they will coexist during the leptogenesis era, 
and there will be lepton flavour violating (but total lepton number conserving) washouts 
processes involving real as well as real A'^i . In what follows we will refer to the washout 
processes involving real or virtual heavy neutrinos as lepton flavour violating washouts 



This implies that, compared to high scale leptogenesis models with Mi > 10^ GeV, 
this electroweak scenario typically suffers from very strong LFVW § . The argument goes as 
follows: the intensity of the LFVW associated with processes involving real or virtual N2 

*In this work we will always consider the case in which the heavy neutrinos are produced thermally 

exclusively through their Yukawa interactions, so that their abundance is null at the beginning of Leptoge- 
nesis. However we have also studied cases in which their initial abundance is that of a relativistic particle 
in equilibrium with the thermal bath (still satisfying Ljv (initial) = 0) and found no significant differences. 

^Similarly strong AL 7^ washouts are also expected in standard (non- resonant) leptogenesis with 
Majorana neutrino masses at the weak scale [12, 13]. 




(2.15) 



(LFVW). 
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is determined by the adimensional parameter rfi2/m* " . For fixed values of the Yukawas 
of N2 - basically for a fixed value of the CP asymmetry in A''i decays- •fh^jra^ scales as 
mpi/M2. Therefore the LFVW will be larger the lower M2. 

Consequently in order to avoid a complete erasure of the asymmetry generated in the 
processes involving iVi, either the CP asymmetry has to be resonantly enhanced so that 
the Yukawas of can be small while still having cnoTigh CP asymmetry, or some of the 
projectors Ka2 have to be very small to reduce the LFVW in some of the lepton flavours. 

Altogether we find the following minimum requirements: 

• The mass of A''i must be not far from the sphaleron freeze-out temperature Mj < 
0{TeV). 

• Either the CP asymmetry is resonantly enhanced or the Yukawa couplings of have 
a strong flavour hierarchy (i.e. Ka2 ^ 1 for a = e, )U or r). 

• Even in the non-resonant case there cannot be a large hierarchy between Mi and M2 
in order to have as little suppression as possible from the factor. 

With these requirements we conclude that to quantitatively determine the viability of 
this scenario we need to consider the evolution of the abundances of both A'^i and (as 
well as the corresponding n'^s) and the lepton flavour asymmetries. In order to do so we 
solve the set of Boltzmann equations which we describe next. 



3. The Boltzmann Equations 

In writing the relevant Boltzmann Equations we first notice that, contrary to the case 
in which the heavy neutrinos are of Majorana nature, the densities of Ni and Ni can be 
different and enter separately in the Boltzmann equations. Thus, in general there is an 
asymmetry between these two degrees of freedom which induces additional washout of 
the lepton asymmetry. Moreover, the usual ALsm = 2 processes mediated by the heavy 
neutrinos (like £ah — > I/jh, etc.) are absent, since total lepton number is perturbatively 
conserved. Therefore, the washout of the lepton asymmetries is due to AL^m = lepton 
flavour violating scatterings mediated by the A^^ {iah — >■ iph, etc.), and ALsm = —ALn = 
±1 reactions with one external iVj. 

Considering all the 1 <-> 2 and 2 <-> 2 processes resulting from the Yukawa interactions 
of the heavy neutrinos and the Yukawa interaction of the top quark, the evolution of the 

^The quantity m, is the equilibrium mass, which is defined by the condition ^^^"'^ = where H 
is the Hubble expansion rate, so that m, = ^^n^''^ ^/gTsM-^ ~ 1,08 x 10~^ eV (g*SM is the number of 
Standard Model relativistic degrees of freedom at temperature T and rupi is the Planck mass). 
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different densities is given by: 

^ lah} + {Niia ^ Qst} + {Nii^ Qs£a} + {NiQs ^ Ua} 

a 

+ J2 {NiNj^iJp} + {Niep^Njeay + {Niea^Njep} (3.1) 

+ ^{NiNj ^ hh} + {Nih ^ Njh}' + {Nih ^ Njh] , 

Y,{Ni ^ £ah) + {Nila ^ Qst) + {Nit ^ Qsia) + {NiQs ^ Ua) 

a 

+ {N.N, ^ Up) + {Niifi ^ NjiaY + {Ni£a ^ N/p) (3.2) 

^ iah) + (Nila ^ Qst) + {Nit ^ Qda) + {NiQ^ ^ tla) 

i 

+ {NiNj ^ iak) + {Niifi ^ NjiJ + {Niia ^ Njip) (3.3) 
+ Yihh^ iJo) + {ifih ^ Uh) + {iph ^ £ah)' , 

where Yx = nx /s is the number density of a single degree of freedom of the particle specie 
X normalized to the entropy density and yx = {Yx — Yx)lY^ (to be used below) is 
the asymmetry density normalized to the equilibrium density. With Qj, and t we denote 
respectively the third generation quark doublet and the top SU{2) singlet. We have also 
defined Y^^ = Yb/3 — Yl^^ , where Yb is the baryon asymmetry and 1/,^ = {2y£^_^ +2/6^^)^'^^ 
is the total lepton asymmetry in the flavour a (with Y'^'^ = Y^^ = Ye^^)- Moreover 
IjVi+iVi = +^jVi is the total normalized density of the heavy neutrino Ni and Yj^._fj. = 
Iat. — Yfj^ is the corresponding Ljv asymmetry. 

To write the Eqs. (3.1)-(3.3) we have also defined the following combinations of reaction 
densities: 



[a, b, . 


■ ■ 1 ear i 
■■ ^ i,J,--\ = -eq-eql \a,b,. 
Ha" 




n j rij . 
■■) eg eoT V') Ji ■•■ 


a,b, ... 


.), (3.4) 


{a,b,.. 


■ ^hj,-) = [a,b,... 


1 - [a,b, .. 






(3.5) 


{a,b,.. 


■ ^ = [a,b,... ^ 


] + [a,6,.. 


■ ^ hj,-] , 




(3.6) 



and the prime written in the contribution of some processes indicates that the on-shell 
contribution to them has to be subtracted. Note that we have not included scatterings 
involving gauge bosons. They do not introduce qualitatively new effects and no further 

density asymmetries are associated to them. We have also ignored finite temperature 
corrections to the particle masses and couplings [26]. In particular we take all equilibrium 
number densities n^, with X Ni,Ni, equal to those of massless particles. 
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-zHs- 



dY, 



Ni+Ni 



dz 



-zHs- 



dY, 



Nj-Nj 

dz 



-zHs 



dz 



After summing over the most relevant contributions " we find: 





-2 


dz 


sHz 




-1 


dz 


sHz 



- 1 



" Ni+Ni 



(3.7) 



7, 



VNi - TFeg Vic + VQz " Vt 

Ni+Ni 



I N^i 



N,Q3 



2yiV, - 2y,„ + I 1 + {VQ, - yt) 

^ ^ Ni+Ni ) 



(3.8) 



dz 



-1 
'sRz 



Ni+N, 



I iVj r 1 I Nil 

i 



p 






Nile 


VNi - 


^ Ni+Ni 
y-eq yia 
Ni+Ni 


Ni+Ni , 


j {yQz - yt) 



+ 




(3.9) 



where we have introduced the notation 7"'^'"' = 7(0, 6, . . . ^ c, d, . . . ). 

In writing Eqs. (3.8) and (3.9) we have accounted for the fact that the CP asymmetry 
in scatterings is equal to the CP asymmetry in decays since only the wave part contributes 
to the CP asymmetry of the different processes (see [27]). Moreover, since the total CP 
asymmetries ej are null, there are no source terms proportional to e^i in the equation for the 
evolution of which is therefore driven only by terms proportional to the different 

density asymmetries. 

It is also important to notice that the equations for the asymmetries are not all inde- 
pendent due to the condition '^^'^-^sm-^n _ where Yb-Lsm-Lm ^ Yr — Ylc., — Yl 



Ylsm = ^^l^^ and Yl^ 



^Yj^._j^.. If we take as initial condition that all the asym- 



metries are null, then Ya„ — Yj^._f}. = 0. 

a i 

In Fig. II we plot the different reaction densities included in the Boltzmann equations, 
normalized to Hn'^'^, where H is the expansion rate of the universe. This normalization is 
appropriate to study the contribution of the different processes to the LFVW. We show 
the figure for Mi = 250 GeV, M2 = 275 GeV, (AtA)ii = 8.2 x IQ-^^ (mi = 10"=^ eV), 
and (A^)22 = 10~^, which are the values of the parameters of one of the examples given 
in the next section. For other values of the Yukawa couplings the reaction densities Jg^Jf^, 



"The scatterings mediated by the Higgs or leptons involving two external heavy neutrinos have been 
neglected because they are much slower than the scatterings involving only one external heavy neutrino 
and the top quark, since the Yukawa couplings of the heavy neutrinos are much smaller than the Yukawa 
of the top quark. 
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10 ^ 1 10 



z = Mi/T 

Figure II: Reaction densities included in the Boltzmann equations as a function of z = Mi/T, 
normalized to Hrip , where H is the expansion rate of the universe. 

7^'^, and -f!t'^^ scale as {X^ while 7^"^ and -ff^- scale as ((AtA)22)^ **• The subtracted 

s-channel reaction density 7^^^ is of the same order and scales as the t- and u-channel 
ones shown in the plot. The infrared divergence of the reaction mediated by the Higgs in 
the t-channel has been regularized using a Higgs mass equal to Mi in the propagator. For 
convenience we have factorized out the flavour projection factors as explicitly displayed in 
the figure. 

From the figure we see that, as explained in Sec. 2.2, the rates of processes involving 
N2 are generically very large (if the CP asymmetry in A'^i decays is required to be not 
too small). When M2 ~ Mi the LFVW induced by processes involving external A''2 are 
dominant over the ones mediated by virtual in the temperature range relevant for 
A'^i leptogenesis. Conversely, if there is some hierarchy between M2 and Mi the LFVW 
mediated by processes involving virtual can become the most important ones, since 
they are not Boltzmann suppressed for T > M2. 

The network of equations (3.7)-(3.9) can be solved after the densities y^^, yh and 
Vt — HQs are expressed in terms of the quantities Ya^, with the help of the equilibrium 
conditions imposed by the fast reactions which hold in the considered temperature regime 

**The contribution of the virtual A'^i to the LFVW processes is negligible in all the cases we are going to 
deal with, so we have only included N2 virtual scatterings. 
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(see [28] and also [25]). For 10*^ GcV '> T > Tc (where we denote by Tc the critical 
temperature of the electroweak phase transition) they read [25] 



m. = -i:</^,^ ^/^ = -E^"^' (3-10) 

13 OL 



with 



221 -16-16\ 
= —(l, 1, 1) and C^ = —\ -16 221 -16 

J 



79' ' 711 



-16 -16 221 



(3.11) 



Moreover, the equilibrium condition for the Yukawa interactions of the top quark implies 
yt - VQs = Vh/'i- 

Above the critical temperature, fast sphaleron processes convert the generated lepton 
asymmetry to baryon asymmetry [4]. Below Tc the Higgs starts to acquire its vev and 
this S'[/(2)/,-breaking suppresses the sphaleron rate, T ^(^b+l) [29-31]. For temperatures 
Mw{T) « T < Mw{T)/aw, T^^b+l) ~ Mw{Mw{T)/awTf{Mw{T)/Tfe^v[-Esp/T] 
[8,30] where aw is the SU{2)l fine structure constant, Mw{T) = gv{^)l\f2 is the W- 
boson mass and the sphaleron energy is Egp ~ M^/i^^^jaw. Because of the exponential 
suppression of ^ i\(^b^l) the lepton asymmetry is not longer converted into baryon asym- 
metry below some temperature Tj for which ^ A{B+L){'^f) / H (Tf) < 1. 

In order to properly account for the evolution of the relevant abundances in the tem- 
perature regime Tf < T < Tc one must extend the system of Boltzmann Equations to 
include the temperature dependent sphaleron rate. The overall effect can be approximated 
by replacing the usual conversion factor ub = ^^'(^b-Lsm) ^ temperature-dependent 
rate given by [31,32] 

77r^ + 54t;(r)^ ^ 
^^^^^ = ^869T2+666.;(r)2^^^«^^)' ^^"^^^ 

where v{T) is the temperature-dependent Higgs vev: 



1 



„(T) = „(l-|)' with r, = „(i + ^ + |i + ^]". (3.13) 

Here A is the quartic Higgs self-coupling, g and g' are the SU{2)l and U{1)y gauge cou- 
plings, and At is the top Yukawa coupling. 

Since the sphaleron processes are effectively switched off at T < Tj, the baryon asym- 
metry is unaffected below this temperature. 

In principle an additional effect is that the set of equilibrium conditions leading to 
Eq. (3.11) are also modified in the temperature range between Tc and Tj. We have verified 
that this effect does not lead to any relevant change in our conclusions. Furthermore for 
T <Tc the effects of the non- vanishing v{T) must be accounted for in the reaction densities. 
As long as Mj is large enough compared to Tc these effects can be safely neglected. 
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4. Results 



In order to quantify the required conditions for generating the observed baryon asymmetry, 
Yb = (8.75 lb 0.23) X 10~^^ [33], we have solved the Boltzmann equations presented in 
the previous section. As discussed at the end of Sec. 2, the CP asymmetry in processes 
involving A'^i or N2 is larger the closer the masses of ^^^^2 and A^i are to each other. Indeed 
the proximity of the masses of the heavy neutrinos is the key parameter that determines 
whether or not successful leptogenesis is possible. We first show how successful leptogenesis 
is possible in this scenario in the resonant mass regime. We then explore the requirements 
for obtaining the observed baryon asymmetry without reaching the resonant condition. 

4.1 Resonant case 

In the case that the CP asymmetry of A^^i decays receives a resonant contribution from A^2 
the proposed mechanism works in a wide range of the parameter space. 
As an illustration we show in Fig. Ill an explicit example in which 



Ml = 800 GeV , M2 = Mi + ^ 



(AU)ii 


= 10-12 


, (AU)22 


= 10-^' 




Kel 


= 0.3, 


K^i = 0.3, 


Krl = 


0.4, 


Ke2 


= 0.1, 


K^2 = 0.1, 


Kr2 = 


0.8 




= pII = 









(4.1) 



With these values the corresponding CP asymmetries are: e^i = 6.5 x 10 ^,6^1 = 3.5 x 
10-2, e^i = -1 X 10-\ ee2 = 1.3 x 10-3,e^2 = 7 x 10"^, and e^2 = -2 x 10"^. 

In deriving the final baryon asymmetry we have taken ruH = 200 GeV for which 
Tc ~ 150 GeV and Tf ~ 100 GeV. 

The figure explicitly shows that Yb-Lsm — as T — > which is mandatory due to the 
conservation of S — L in this scenario and the assumed initial conditions. Indeed it can 
be verified that, although Y/^^ saturate at a finite value at low temperature, with the sign 
assignments for the asymmetries, at any z it is verified that Yb-Lsm ~ l-^Ael + I^A^I ~ 
\YaA = Yif, = \Yj^^_j^\ + ^ |5jv2_]^|. Still as illustrated in the figure, the 

observed baryon asymmetry is generated once the sphalerons switch off below Tf = 100 
GeV {z>8). 

We notice that despite the CP asymmetry e^i » 60,2, the asymmetries verify 
l^^-Tv^l ^ I^A'2-JV2l' "^^^^ ^'^ because, even though the lepton asymmetries y^^ are 
mostly produced in processes involving Ni, it is the inverse decay £ah Ni what deter- 
mines how much of the lepton asymmetries is transferred to the iVj asymmetry. Moreover, 
recall that there are no source terms proportional to eai in the evolution equations of 
Yj^__j;f. Eq. (3.8). Then, since the N2 Yukawa couplings are larger, the inverse decays of 
the A''2 are more efficient and a larger A''2 asymmetry is produced. Consequently for the 
small A'^i Yukawa couplings considered, the A'^i asymmetry plays a very little role in the 
dynamics of the system. 
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Figure III: The asymmetries \Ya^ \ , \ Y]^._^. \ , Yb-Lsm j and tlie densities Yjq.j^jq. , Y^^^ as a fmic- 
tion of z for the values of the parameters given in Eq. (4.1). 

We note also that contrary to leptogenesis scenarios which occur well above the elec- 
troweak scale, once the heavy neutrinos have decayed, the universe is left with an equal 
amount of lepton and baryon asymmetry. Furthermore the flavour asymmetries Yao typi- 
cally remain some orders of magnitude greater than Yb- 

4.2 Non-resonant case 

In this regime, in order to have enough CP asymmetry at least one of the Yukawas of 
has to be very large. The strongest bounds on the Yukawas of the heavy neutrinos 

for the range of masses we are interested in come from constraints on violation of weak 

universality, lepton flavour violating processes and collider signatures [21] which allow for 

|(AtA)22| <5x 10-3(M2/f)2. 

However for flavour effects to be relevant, the Yukawas of the N2 must be smaller than 

the Yukawas of (at least) one of the charged leptons, so we require that X2a ^ A,- ~ 10~^ . 

This is so because if the fastest leptonic reaction rates were those associated to the 

^^This is a rough estimate of the vahdity of the "fuUy flavor regime" which is enough for our purposes; 
for a more detailed analysis about this point see [34]. 
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Figure IV: The asymmetries I^a^ |, I^Wi-iVi L ^s-isM j the densities ^/Vi+JVi: ^jv'+jv ^ fmic- 
tion of z for the values of the parameters given in Eq. (4.2). Note that the densities Y^.^^. and 
^t^ .-fr. have been rescaled by 10~^ in order to fit into the figure. 



flavour basis which diagonalizes the density matrix would be that formed by the lepton in 
which decays, ^2 = ^a'i^a/ X^q, I'^a2p) together with two states orthogonal to £2 and 
hence no baryon asymmetry would be generated since in this basis the CP asymmetry is 
zero. 

Furthermore if the Yukawa couplings of N2 with all the three light leptons were com- 
parable and large (even if smaller than At) there would be strong LFVW in all flavours 
(see Fig. II) and the generation of the baryon asymmetry would be strongly suppressed. 
Therefore, as was explained before, some of the Yukawa couplings have to be very small. 
Note that for a fixed value of (A^A)22 the CP asymmetry in a given flavour decreases lin- 
early with decreasing Yukawa coupling of that flavour with A'^2! while the washout terms 
decrease quadratically with it. We conclude that the most favorable situation is to have a 
strong hierarchy in the flavour structure of the Yukawas. 

Concerning the optimum range of Yukawa couplings for A^i there are two relevant 
effects. On one hand if they are too large - fhi > ~ 10~^ eV - the maximum values of 
the B — LsM asymmetry occur at z < 1, so that Mi < Tj in order for the sphalerons to 
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decouple when the baryon asymmetry is largest. In this case the analysis is more complex 
because one cannot neglect the effects of the breaking of SU (2) in the reaction densities. 
The expected effect is the reduction of the Ni decay rate due to the phase space factors. 
In order to determine the effect on the final B asymmetry a dedicated study is required 
which is beyond the scope of this paper. On the other hand, for values m\ ^ m* , the peak 
in the B — Lsm asymmetry shifts to large values of z but it is in general lower because 
of the smaller production of A^i when starting from a zero abundance as initial condition. 
One may wonder if this conclusion may be modified when assuming a non- vanishing initial 
Ni abundance which would allow for very late A''! decay. It is not, because at large values 
of z the LFVW are very suppressed and therefore the flavour effects, which are essential 
in this scenario, do not survive. In summary, generically larger baryon asymmetries are 
expected for mi ~ m*. 

With the above considerations in mind, we have explored the parameter space for a 
fixed value of M2/M1 (chosen near to 1) and in Fig. IV we plot the evolution of the different 
asymmetries and densities for a set of parameters representative of the cases with highest 
production of baryon asymmetry: 

Ml = 250 GeV, M2 = 275 GeV , 
(AtA)ii = 8.2 X 10-^^ (mi = 10"^ eV) , 
(AtA)22 = 10-^ , 

^ ' ' (4.2) 

Kei = 0., K^i=0.3, Kri=0.7, 

Ke2 = 0. , Kf,2 = 10-^° , Kr2C±l, 

12 _ 12 _ 12 _ 1 
PefM Per Pfir • 



Since ^/{)'JXj22 > we have safely chosen the projectors K(,i,Ke2 = to prevent any 
possible flavour projection effects associated to the Yukawa interactions of N2, which would 
complicate the description of the problem without substantially changing the results. In an 
eflFectively two flavour case the CP asymmetries are proportional to only one phase factor, 
in this case to p^^, therefore we have adopted in the example its maximum possible value. 
With these values the corresponding CP asymmetries are: Cgi = ee2 = 0, e^i = — e^-i = 
8.7 X 10-i\ e^2 = -er2 = 8.7 x IQ-^i. 

Prom the flgure it can be seen that even with these large N2 Yukawa couplings, their 
strong flavour hierarchy and the small l/(a2 — 1) suppression, the produced Yb falls short 
to explain the observations by about 5 orders of magnitude. However, we notice that in this 
regime the asymmetry comes mainly from processes involving Ni . Therefore the B — Lsm 
asymmetry is approximately proportional to l/(a2 — 1) which is around 5 in the example. If 
we take M2 closer to Mi , that factor and the corresponding Yb grow accordingly. Thus we 
see that in this regime it is also possible to generate the required baryon asymmetry as long 
as Ni and N2 are strongly degenerated even if still not in the resonant regime. For example 
for the values of parameters given in Eq. (4.2) the observed baryon asymmetry could be 
produced if 02 — 1 ~ 2.4 x 10~^, which is still far from the resonance ((72 = 4 x 10~^). 

Something to note is that despite the great hierarchy among the projectors of onto 
£^ and ir, the flavour asymmetries Yac* (« = t) quite similar in size. This is because 
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the evolutions of the different asymmetries are strongly coupled due to the conservation of 
B — L and the null value of the total CP asymmetry in A^i decays. 

To obtain an estimate of the order of magnitude of the density asymmetries as well as to 
understand their dependence on the N2 flavour projectors (which are more relevant than the 
Ni projectors when (A'^A)22 S> (A^)ii), we have developed the following semiquantitative 
approximation: 

(i) From Fig. II we see that for K^2 ^ 10~^ the rates of the processes N2 ^^.r/i are 
larger than the expansion rate of the Universe in the most relevant range of temperatures. 
Therefore at each instant the thermal bath has time to relax, i.e., the production of asym- 
metry equals its erasure. This implies that the derivatives of the density asymmetries are 
negligible with respect to the source and LFVW terms, hence we set them to zero in the 
Boltzmann equations. 

(ii) We keep the CP asymmetries produced by A^i (since £0,2 ^ f-ai), as well as the N2 
decay and inverse decay reactions, and neglect all the remaining subdominant contributions, 
including the partial conversion of lepton asymmetry into baryon asymmetry during the 
leptogenesis era. 

Within this approximation, Eqs. (3.8) and (3.9) simplify to 

S,{z) + K^27^^'[yr,,-ye,] =0, (4.3) 
5.(z) + K,27ilMyiV.-2/d =0, (4.4) 

where Sn{z) = —St-{z) = e^i ( Jlj"*"^^ — l) 27^^^ is the source term normalized to {—sHz)~^ 

A third equation is provided by total lepton number conservation, i.e. Y/V2-JV2 + 
+ = 0. In the most relevant temperature range for the model, T ~ M2, we can 

approximate the equilibrium density of N2 by that of one relativistic degree of freedom, 

therefore 

VN, + 2ye, + 2ye, = . (4.5) 
When K^2 *C 1 ^ Kt2 the solution to this system of equations is 

yN2 = Vir 

From this analysis, it is clear that despite the large hierarchy in the projectors of 
A^2, all the density asymmetries have the same order of magnitude. Moreover, since the 
source term is proportional to y^^^, the density asymmetries are inversely proportional 
to ^JW^. We have verified numerically that this dependency of the asymmetries on the 
projector actually holds in the range 10~^ < K^ja < 10~^, where the approximations we 
made are expected to be valid. 

For ^ K^2 ^ lO"^'^ the rates of the processes N2 i^h are lower than the expansion 
rate of the Universe, hence point (i) is not longer true. In this range of K^2, the density 



3 S^jz) 
2 S^{z) 



N2 



(4.6) 
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asymmetries decrease as \fK^ because the main dependence on i^^2 comes from the 
relation e^^ri oc ^J~K^. Thus, fixing all the parameters but K^2 to the values given in 
Eq. (4.2), the baryon asymmetry is maximized for K^^i between 10~^ and 10~-^°. 

5. Summary 

In this article we have studied the possibility of generating the observed baryon asymmetry 
via leptogenesis in the decay (and scatterings) of heavy Dirac Standard Model singlets with 
O (TeV) masses in a framework with B — L conservation above the electroweak scale. In 
this scenario a total lepton number, which is perturbatively conserved, can be defined. 
This lepton number is shared between the Standard Model leptons and the heavy Dirac 
singlets as described in Sec. 2. In this scenario, despite the total CP asymmetry is null, a 
CP asymmetry in the different SM lepton flavours can be generated (sec Sec.2.1). 

The additional physical condition for generating a non- vanishing B in this framework 
is that the sphalerons depart from equilibrium during the decay epoch. For symmetric 
initial conditions (no net baryon nor total lepton number present), B — LgM = after the 
heavy neutrinos have disappeared. Consequently if the sphalerons were still active after 
the heavy Dirac singlets decay epoch, the final baryon asymmetry, being proportional to 
B — LsM, would be zero. However if they depart from equilibrium during the decay epoch 
the baryon asymmetry freezes at a value which in general is not null. 

In summary in this scenario the baryon asymmetry is generated by the interplay of 
lepton flavour effects and the sphaleron decoupling in the decay epoch. In order to quan- 
tify whether enough baryon asymmetry can be generated we have constructed and solved 
the network of relevant Boltzmann Equations associated with the abundances of the two 
lightest Dirac heavy singlets, their asymmetries and the three SM flavour asymmetries. 
The results are given in Sec.4. 

We find that the ratio of the masses of the two heavy Dirac neutrinos is the key pa- 
rameter that determines whether or not successful leptogenesis is possible. The relevant 
Yukawa couplings are constrained from above by the requirement of having flavour ef- 
fects, and from below by the requirement of large enough CP asymmetry. Within these 
boundaries we conclude that successful leptogenesis can occur if the masses of two heavy 
Dirac singlets are very degenerate M2/M1 — 1 < 0(10^^). Recall that in our framework 
this degeneration is not a consequence of the small breaking of total lepton number, as 
in other low scale leptogenesis scenarios [8,10], although it may be due to an additional 
symmetry of the heavy sector. In particular, if the CP asymmetry is resonantly enhanced 
- (M| — Alf) A!iT2 - , the mechanism works for a wide range of values of the Yukawa 
couplings and flavour projections. It is worth to explore whether the heavy neutrinos will 
be observable at the LHC and/or lead to measurable lepton flavour violating signals, within 
the parameter regions allowed by leptogenesis. 
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